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Abstract 

Properties related to entanglement in quantum systems, are known 
to be associated with distinct properties of the corresponding classical 
systems, as for example stability, integrability and chaos. This means 
that the detailed topology, both local and global, of the classical phase 
space may reveal, or influence, the entangling power of the quantum 
system. As it has been shown in the literature, the bifurcation points, 
in autonomous dynamical systems, play a crucial role for the onset of 
entanglement. Similarly, the existence of scars among the quantum 
states seems to be a factor in the dynamics of entanglement. Here we 
study these issues for a non-autonomous system, the Quantum Kicked 
Top, as a collective model of a multi-qubit system. Using the bifurca- 
tion diagram of the corresponding classical limit (the Classical Kicked 
Top), we analyzed the pair- wise and the bi-partite entanglement of 
the qubits and their relation to scars, as a function of the critical pa- 
rameter of the system. We found that the pair-wise entanglement and 
pair-wise negativity show a strong maximum precisely at the bifurca- 
tion points, while the bi-partite entanglement changes slope at these 
points. We have also investigated the connection between entangle- 
ment and the fixed points on the branch of the bifurcation diagram 
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between the two first bifurcation points and we found that the entan- 
glement measures take their extreme values precisely on these points. 
We conjecture that our results on this behavior of entanglement is 
generic for many quantum systems with a nonlinear classical ana- 
logue. 

PACS : 03.67Mn, 05.45Ac, 05.45Pq, 05.45Mt 

KEY WORDS : Entanglement, Quantum Chaos, Bifurcations, Scars, 
Kicked-Top Model. 
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1. INTRODUCTION 



Quantum systems with classical analogues can show properties which 
strongly depend on the detailed topological structure of the classical phase 
space. One of the oldest and intensively studied properties is the quantum 
signature of classical chaos, the so called quantum chaos |T]. In a broad 
sense , the last concept is mostly defined in relation to the classical system. 
Namely, knowing the parameter ranges for which the corresponding classical 
system is chaotic, one associates quantum structures with characteristic clas- 
sical behavior in these parameter ranges, as for example the well established 
eigenvalue and eigenvector distributions [3]. On the other hand it has 
been a lot of research activity for a purely quantum mechanical definition 
of chaos in terms of intrinsic quantum properties, the most characteristic 
of these being the quantum entanglement [1]. Indeed it has been shown 
that the entangling power of the evolution operator and the stability of en- 
tanglement are strongly influenced by the properties of the classical regime 
Pia[l0l[IIl[I3[7l[8l|13l|lll|15lll6lll71. Though there are some conflict- 
ing results and mostly dependent on specific models, the consensus is that 
quantum entanglement is an effective signature of quantum chaos. One then 
may pose the question whether the more detailed structure of the classical 
phase space may reveal a corresponding structure or behavior of the quan- 
tum system. Thus it is natural to be asked if the bifurcation points can 
signal anything for the quantum regime. After all, the most common scenar- 
ios of classical chaos are associated with specific series of bifurcations. The 
question that has been analyzed in the literature is the behavior of entan- 
glement at the bifurcation points, and it has been shown ,for autonomous 
systems, that indeed entanglement is sensitive to this classical information 
[T8l [T9l [20} [2T] . In a different direction, there is an established relationship 
between classical periodic orbits and localization of wavefunctions, the so 
called scars [221 1231 [211 [25], [26l [27] . The existence of scars is closely related to 
the chaotic regime, and an immediate question is whether entanglement and 
scars are related . There exist studies of autonomous systems which answer 
this question positively [31]. In this paper we study these issues for a non 
autonomous system, the Quantum Kicked Top [28l [291 [20]. We use this to 
model the collective behavior of a system of qubits perturbed externally by 
discrete kicks. This system has been studied extensively in the literature. It 
is well understood how the classical and quantum chaos are interrelated, what 
is the role of periodic orbits, how entanglement production may be quantified 
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and analyzed with this model and finally how entanglement may be used as 
a signature of quantum chaos. Here first we study the role of bifurcations 
in the entangling power of the kicked top. We have found that the pair-wise 
entanglement shows a strong maximum precisely at the bifurcation points. 
Then we analyzed the relation of the existence of scars to entanglement prop- 
erties and we have observed that the bi-partite entanglement as a function 
of the critical parameter shows a strong change of slope. Finally to have a 
more complete picture, we investigated the dependence of entanglement on 
the position of the phase space point with respect to the branch of stable 
fixed points. We found that , precisely on this branch , between the first two 
bifurcation points, entanglement becomes extremum. Of course a complete 
understanding of this behavior would come from an extensive analysis of the 
relation between entanglement and the stability properties of the points on 
higher branches of the bifurcation diagram, high values of k, well inside the 
chaotic regime. In section 2 we introduce the model and review some known 
facts from the literature. We also present the concept of quantum scars as 
it is applied in this model. In section 3 we present our calculations and in 
the final section we discuss our results and present a conjecture about their 
generality. 

2. BACKGROUND 

2.1 Quantum and Classical Kicked Top 

For a system of N qubits we define the collective operators 

Ja = T.-y (1) 

i=l ^ 

with {(Tja} the Pauli matrices. We use as a basis the symmetric Dicke states 
[33] m >,m = —j, ■ ■ ■ ,j,j = N/2}. This means that the dimension of the 
space is N+1 and not 2^, a fact that facilitates the numerical computations. 
We model the dynamics of this system with the Quantum Kicked Top (QKT) 
Hamiltonian [23j. This is the non autonomous quantum system 

+00 

H{t) = {hp/T)Jy + {hk/2j)Jl Sit-nr) (2) 

n=— oo 

The linear term describes a precession around the y axis with angular fre- 
quency p/t. The quadratic perturbation , with perturbation parameter k 
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accounts for periodic kicks in time steps with period r. The unitary time 
evolution is given by the Floquet operator 

with the normahzation h = 1,t = 1. The evolution is a precession around 
the y axis by an angle p with a subsequent kick. If the initial state of the 
system is \yo >, after n units of time the state is given by |y„ > = -F"|yo >■ 
The equation 

F|$^>= e-'"|$^> (4) 

determines the pseudo-eigenfunctions and the pseudo-eigenvalues of the Flo- 
quet operator. These functions provide a complete set for the expansion of 
the states of the system, namely 

|yn $^e^"^"'|$^ >< ^m\yo > (5) 

m 

We note here that this dynamics leaves the set of symmetric states invariant. 
The transition to the corresponding Classical Kicked Top is achieved by 
first expressing the dynamics of the system in the Heisenberg Picture. The 
evolution of the angular momentum operators 

4 = Ft^JoF" (6) 

gives the equations 

i; = -iXcosp + J.sinp + J^)e^f{-^^co.p-j..mp+i) ^ 

j'. = hiJxCosp + J.sinp + iJJe*j(-^-'=°-f-^---f+5) + h.c (8) 

Jj. = JzCosp — JxSinp (9) 

In order to define the corresponding Classical Kicked Top one has to take 
the limit of very large angular momentum j. Defining the new operators 
fa = Ja/j , a e {x, y, z} which satisfy [fa, ft] = '-^fc, we have that for 
J — oo these operators commute. Then writing X = T^, Y — Ty, Z — Tz, in 
this limit , we have the classical equations of motion 

X' = Re{Xcosp + Xsinp + iY)e'''^^^°'P-^'''^P^ (10) 
Y' = Im{Xcosp + Xsinp + iY)e'''^^''°'P-^''''P^ (11) 
Z' — —Xsinp + Zcosp (12) 
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For the specific value p = 7r/2 the dynamics is given by the equations 



X' 
Y' 
Z' 



Zcos{kX) + Ysin{kX) 
-Zsin{kX) +Ycos{kX) 
-X 



(13) 
(14) 
(15) 



and since we have 




X^ + Y^ + Z^ = 1 



(16) 



the transformation X = sinOcoscj) , Y = sinOsincj) , Z = cos6, gives a two 
dimensional map as the equations of motion for the Classical Kicked Top. 
This classical dynamics is well analyzed in the literature [281 EHl EQ]. For 
small values of k we have islands of regularity which are destroyed as k in- 
creases. For intermediate values of k we have a mixed phase space, while for 
big values we have full chaos. There are two trivial equilibrium points, at 
the north and south poles of the sphere, defining these points here w.r.t. Y 
axis. We follow the bifurcation diagram of one of these equilibrium points 
analyzing the behavior of the corresponding quantum system at the bifurca- 
tion point. In a more complete analysis , we plan to track the signatures of 
many bifurcations points and picture the overall significance of the bifurca- 
tion landscape for the quantum behavior. There are generic scenarios of the 
onset of chaos, and these are well studied and understood, as well as cases 
of systems with singular behavior. We consider the study of the quantum 
signature of these classical portraits as very important. In this paper we are 
only taking the very first step of such an analysis. 

2.2 Quantum Scars. 

In the semi-classical approximation theory of Quantum Mechanics the 
classical periodic orbits play a fundamental role |lj. It has been proved 
that the fluctuating part of the density of states can be written as a sum 
of terms over the classical periodic orbits, and that the existence of these 
orbits influences the energy spectrum of the corresponding quantum system. 
But one of the most impressive and unexpected results is the existence of 
scars [22|. These are states which are strongly localized in the phase space 
at points which are directly related to the classical periodic orbits. These 
Quantum Scars (QS) do not exist for every periodic orbit and not every 
quantum state is related to any periodic orbit. But it is interesting to inves- 
tigate how the existence of QS affects the quantum properties of the system. 
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For non-autonomous systems, like the QKT, scars can be found among the 
pseudo-eigenstates of the Floquet operator. It has been proved that scars 
play a significant role in the time evolution of a coherent state which starts 
from a point of the periodic orbit |22l |23l [2ll |25l [261 1221 • In this paper we 
analyze the relation of scars to the entanglement. We base our analysis on 
some known results. 

Let I7 >= \9, > be a given state of the QKT where \9, > is the spin 
coherent state defined by i32|, {|^, >= R{9, >; -vr < < vr, < < 

tt}, with 

R{9,(j)) = exp{i9[JxSin(f) — JyCoscj)]} (17) 

Then for the nth iterate of the Floquet operator we have I7' >= F^\'-f >. 

Here 7' is the coherent state parameter given by the classical dynamics 7' = 

T(- ■ ■ T(T(7))), where T is the appropriate classical limit of the Floquet map. 
^ ^ 

n 

Thus for a classical periodic orbit of order K we have 70 = T(- ■ ■T(T(7o))), 

V ' 

K 

[28] . Then for the M-step propagator /m =< 7o|-F*^|7o > we have, after 
expanding in the pseudo-eigenfucntion basis 

N 

fM = E^n(7o)e^'"^" (18) 

n=l 

where N is the dimensionality of the system and if„ = | < $n|7o > P 
the Husimi representation [23] of the nth pseudo-eigenf unction The 
importance of this quantity is that its Fourier Transform 

L N L 

L= y: ^-'^''fM = E ^"(70) E e'""'^'"""^ (19) 

M=-L n=l M=-L 

contains the information for the quantum scars. That is, when u is close 
to a pseudo-eigenenergy, the spectral density is dominated by the overlap 
function of the corresponding pseudo-eigenfunction if the latter is associated 
to a scar. Thus we locate the QS that corresponds to the periodic orbit 
r = (70) 7i) ■ ■ ■ ? 7Af-i) as a maximum in Stated differently, the eigen- 
function associated with the energy at the peak is scarred by this periodic 
orbit [23]. We refer to the literature for more details on scars and their prop- 
erties. 

2.3 Entangling Measures. 
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Our main objective is to study the pair-wise and the bi-partite entagle- 
ment in a system of N qubits which is subjected to a time-dependent external 
perturbation and analyze the interrelation of these quantities with the ex- 
istence of bifurcation points of the corresponding classical system and the 
role of quantum scars. As we said above, we model this collective behavior 
as Classical and Quantum Kicked Tops, the properties of which are very ex- 
tensively analyzed in the literature. We quantify the pair-wise entanglement 
with both the entanglement of formation and negativity and the bi-partite 
entanglement with the linear entropy. More specifically, let p be the density 
matrix of the system and \i1i2 ■ ■ - in >ik € {0, 1} the computational basis. 
Let pk = trALL^kip) and pkm = ALL^k,m{.p) be the one qubit and two qubits 
reduced density matrices. We use three measures of entanglement. The bi- 
partite entanglement Q = 2 — -|: X]i tr^pf), the Negativity Np = N^pk^n) and 
the Entropy of Formation Ep = EformationiPkm)- Negativity is defined as 
follows : For a generally mixed state p, N{p) = J2j fnax{0, —pj} where {pj} 
are the eigenvalues of the partial transpose of p. Entanglement of Formation 
is defined as Ef = H( ^^^^~ —) where H(x) is the Shannon entropy function 
and C the concurrence defined by C = max{0, Ai — A2 — A3 — A4} and the 
quantities Xj are , in decreasing order , the square roots of the eigenvalues of 
the matrix pi2{(^iy ® a2y)p*i2(^ij/ ® (^2y))- We refer to [Ml EZl ESI SD ESI SD] 
for more details on the properties of these quantities. In the Discussion we 
comment on the use of Generalized Entanglement [3^ [35| [T7] . 

3. THE RESULTS 

3.1 The Classical Picture. 

We first present a series of portraits of the Classical Kicked Top for different 
values of the parameter k and the bifurcation diagrams for the angles 60 and 
00 of the equilibrium point of the south pole. 

In Figure 1 we reproduce the phase space portrait for four different values 
of the coupling constant parameter k, with one value at the bifurcation point 
k=2. We present these four figures to point to the way the classical phase 
space picture is evolving with the increase of parameter k. Since we are not 
studying here the classical behavior of this system, we are not concerned with 
the small structures evident in the figures as we approach the chaotic regime. 
It is , of course, very interesting to study the relation of these details to the 
quantal properties of the associated system. We only note here, that in the 
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limit of fully developed chaos, there exist in general unstable periodic orbits 
which we expect to play a role in the behavior of the quantum system. In 
Figure 2 we plot the bifurcation diagrams for the angles 9 and which start 
from the equilibrium point in the south pole. We note that the fuzziness of 
the bifurcation diagrams after the point k=2 and before the area of complete 
chaos is due to numerical inaccuracies. In this figure, the first bifurcation 
point is more than evident. In our results , below, for the entanglement, 
the signature of the existence of this point in the behavior of the quantum 
system is clearly seen. But it is not yet clear to us how to pin point the onset 
of the band. Perhaps, as it is commented in the Discussion, other measures, 
like Generelaized Entanglement may offer a better picture. 

3.2 Bifurcation Points and Entanglement of Scars. 

As it was described in the Introduction, scars are quantum states which 
are associated to classical periodic orbits and appear as prominent structures 
in phase space through their Husimi representation. In order to relate the 
entanglement of the scars to the bifurcation point at k=2, we first identify 
the scar located at the south pole for each value of k, from k=0 to k=6.8. 
Then for each value of k we compute for the corresponding scar the pairwise 
negativity, the pairwise entanglement of formation and the bipartite entan- 
glement. We did the calculations for j=35, that is N=70 qubits. In Figure 
3 we plot the scars and the associated spectrum corresponding to the values 
of k in Figure 1. In this figure it is seeing that as k increases more peaks 
arc present, and as k approaches the region of full chaos the peaks multiply 
greatly. We identify the highest peak with the scar that we study. The cor- 
responding plot of the Husimi function shows a progressive delocalization. 
In Figure 4 we plot the three different entanglement measures, pair-wise 
negativity, pairwise entanglement and bipartite entanglement for N=50 and 
N=110 qubits. First we observe that in all plots the signature of the bifur- 
cation point at k=2 is more than evident. Second, there is not any dramatic 
change with the increase of the number of qubits. In this figure we see that 
though the entanglement is associated clearly with the scarred function at 
the bifurcation point, this association is not strict, and it shows a relatively 
wide span of values where the entanglement is appreciable. It seems that 
this is related to the progressive delocalization of the scarred state. 

3.3 Entanglement and Periodic Orbits. 

In Figure 5 we plot the pair- wise entanglement of formation and the bi-partite 
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entanglement for time evolved states. Specifically, we start with two coher- 
ent states located at the south and north poles respectively, and we compute 
their entanglement after 500 kicks. We repeat this for a large number of k 
values from k=0 to k=6.8. We observe again the trace of the classical bifur- 
cation point at k=2. There is not any difference between south and north 
poles and not a strong dependence on the number of qubits for N in the 
order of 10. It is interesting to note that classically the two poles behave 
differently at the bifurcation points |28|. Here , again, though the value of 
k=2 is evidently seeing to be prominent, there is not a sharp peak. It seems 
that the bifurcation point is smoothly felt by the periodic orbit starting from 
the equilibrium points. A more extensive study is probably needed to under- 
stand the influence of the local structure of the classical phase space on the 
quantum behavior. 

3.4 Entanglement along the branch of stable points. 

To have a more complete picture of the role of the classical stable fixed 
points we analyzed the dependence of entanglement on the location of the 
fixed points w.r.t. the branch of the stable points between the values k=2 
and k=12.73, the latter corresponding to the second bifurcation point. We 
computed the bipartite entanglement and the pairwise negativity for vari- 
ous values of k in the above interval. For the k=12.73 the fixed point is 
{9, (f)) = (2.32, 0.38), and we computed the dependence of entanglement on 9 
and (f). We found that, at precisely this fixed point, the entanglemnt is min- 
imum. In Figure 6 third row, we plot the bipartite entanglement vs either 9 
or (j). In the same figure, first and second rows we plot the results for k=2.2 
and k=2.4. For these cases the corresponding values for the fixed points are 
{9, 0) = (1.9, 1.2), (2.1, 1.0). We keep 9 fixed and we plot both the bipartite 
entanglement and the pairwise negativity vs 0. Again we observe that the 
fixed points give minimum entanglement. We note that the same behaviour 
is observed for all values of k in the interval [2,12.73]. 

4. DISCUSSION 

In this paper we posed the question of the relation of the local struc- 
tures of classical phase space to the purely quantum mechanical property of 
entanglement. Generally, it is well established that, both qualitatively and 
quantitatively, entanglement and its fragility, namely the properties of the 
decoherence process, strongly depend on the classical regime. That is, the 
presence of classical structures, in the corresponding classical system, can be 
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traced in the properties of the quantum system. But the onset of classical 
chaos is generically related to biftuxations from equilibrium points, and the 
scenarios of chaos are pictured in connection to local or global changes of 
the phase space portrait of the system. It is generally expected that a global 
change of the classical landscape can be felt by the quantum system. But it 
seems very natural to ask whether a local property, that is a single bifurca- 
tion point can be traced in the quantum properties. It is this question that 
we have tried to answer here using scars and periodic orbits as "tools" of 
the analysis. We have used the Classical and Quantum Kicked Tops because 
they have been well studied and understood. We use the top as a collective 
model of N qubits. We applied three measures of entanglement , pair-wise 
negativity, pair- wise entanglement of formation and bi-partite entanglement. 
We traced the behavior of the system for many values of the coupling pa- 
rameter k, from k=0 to k=6.8. As it is known for k=2 we have the first 
bifurcation point of the equilibrioum at the south and north poles. First we 
presented four portraits of the phase space for characteristic k values and the 
bifurcation diagrams of the south pole. Then for each value of k in the above 
range we locate the scarring state and compute the entanglement. The scars 
along with the associated spectra are plotted for the selected four values of 
k. Finally we propagated two coherent states , initially located at the south 
and north poles and we computed their entanglement after 500 kicks. In all 
plots of the measures of entaglement the bifurcation point at k=2 is more 
than evident. We conjecture that this is not a non-generic behavior. That is 
we expect that the bifurcation points are intimately related to the quantum 
structures. But , as it happens in the case of scars, we cannot tell whether all 
bifurcation points play this role and moreover whether the sensitivity of en- 
tanglement can always be related to a bifurcation phenomenon in the classical 
regime. To establish a more complete picture, at this stage of our investiga- 
tion, we made an analysis of the dependence of entanglement on the propery 
of stabihty of the fixed points. On the branch of stable fixed points in the 
bifurcation diagram between k—2 and k=12.73 (first and second bifurcation 
points), for various values of k, we computed the dependence of entanglement 
on the position of the point in phase space. We observed that when the point 
crosses the branch the entanglement becomes minimum. We ploted the cases 
for k=2.2, 2.4 and 12.73. The behavior is similar for other values of k in this 
interval. Given the fact that the bifurcation diagram is very complicated 
it would be interesting to see whether for very high values of k, we have a 
new pattern on the entanglement dependence. We are currently investigating 
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these issues. We end with few comments on the use of Generahzed Entan- 
glement (GE) [511 ESI [IZ]- The static or dynamic signatures of quantum 
chaos, hke eigenvalue distributions or fidelity, and entanglement generation, 
with the associate entanglement power of the evolution operator are widely 
used, but depend strongly on the basis employed. But more importantly, 
they depend on the partition of the given system into subsystems, or the em- 
bedding into a larger system. It is this practice which we follow here using 
the pair-wise and the bi-partite entanglent and the negativity. On the other 
hand the GE depends on a selected subalgebra of observables. The regime 
is characterized by sums of averages of these operators. These averages are 
sensitive detectors of the global classical phase space, and for this reason GE 
is considered presently as the most effective signature of quantum chaos. In 
the present paper we were interested on two critical details of the classical 
dynamics, namely the role of bifurcation points and of general periodic solu- 
tions in the entanglement affected by the dynamics on the subystems. Thus 
from the start our interest is on subsystems and their combined state, and 
according to our bi-partite entanglement measure is actually related to 
GE. Of course we consider the question of the dependence of GE on the local 
properties of the phase space as very important, and we are looking into this 
issue. 
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Figure 3l ScarE for the south pole equiUbriuin poiut. First row k=1.5jState 
n=4T with pseudoenergy ijjVi=4.2922, Second row k=?, 11=47. u^n =4,05031, 
Third row k-2.5,n=51 ;4j„=4.3933. Fburth row k=4.5, n=45 ij,i= 3.8732 





J@?teii^)eiHiErit measures for N=[50 quIbrfe-^iiaoBth pole (left) m^' 
jf=il&.i9^ib&B-Q.o(rth pQle(right). Firsli row pairwLse .iiegtdii«^f:--:&KMiEd-T?^ 




Figure S: EkitaiLgfemeEl; of the orbitB tor N=?3t) qilbits, startmg &0111 tlie.sosltli' 
;:pole (left) and N=70 qubits. starting &om the north pofo: (i^^Jii- i^rafeKSW. 



6=1. e p4S^ 





f^ul^^: Entanglement along tbe braiich of stable points. Filst row bipatiitft.. 
s;'nLui3;k'iiK']it tiiid p.dr'hvisi? ]if'i3;iiti\'ity fui: k=2.2. Sc-i/'uinl lov' i.lit? s^ap^-jSsff. 
k=3 ,4, Third row bij>artite antangleuieut for )i= 12,73, tte v&lue for .tbeipeEOEd' 
bifurcation. t* 



